Asymptotically Friedmann self-similar scalar field solutions with potential 

^Masanori Kyo, ^Tomohiro Harada* and ^■•^Hidcki Maeda^ 

^Department of Physics, Rikkyo University, Toshima, Tokyo 171-8501, Japan 

^Centro de Estudios Cientificos (CECS), Arturo Prat 514, Valdivia, Chile 

^Department of Physics, International Christian University, 

3-10-2 Osawa, Mitaka-shi, Tokyo 181-8585, Japan 

(Dated: June 16, 2008) 

We investigate self-similar solutions which are asymptotic to the Friedmann universe at spatial 
infinity and contain a scalar field with potential. The potential is required to be exponential by self- 
similarity. It is found that there are two distinct one-parameter families of asymptotic solutions, 
one is asymptotic to the proper Friedmann universe, while the other is asymptotic to the quasi- 
Friedmann universe, i.e., the Friedmann universe with anomalous solid angle. The asymptotically 
proper Friedmann solution is possible only if the universe is accelerated or the potential is negative. 
If the potential is positive, the density perturbation in the asymptotically proper Friedmann solution 
rapidly falls off at spatial infinity, while the mass perturbation is compensated. In the asymptotically 
quasi- Friedmann solution, the density perturbation falls off only in proportion to the inverse square 
of the areal radius and the relative mass perturbation approaches a nonzero constant at spatial 
^ ' infinity. The present result shows that a necessary condition holds in order that a self-gravitating 

^ ' body grows self-similarly due to the constant accretion of quintessence in an accelerating universe. 
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I. INTRODUCTION 



It is now widely believed that the expansion of our universe got accelerated in its early phase of evolution, which 
bi)! is called inflation. Under reasonable assumptions, this inflation implies that some form of matter fields with largely 
negative pressure may have dominated the energy of the universe. The simplest inflation models postulate that a 
-^ I scalar flcld with flat potential would have induced this early-phase acceleration. On the other hand, the independent 
^ . observations of supcrnovae, the cosmic microwave background, and the large scale structure has recently revealed 
f^ ■ that our universe is currently in a phase of accelerated expansion [l|, Q • This implies that the energy of our universe 
is currently dominated by some form of matter fields with largely negative pressure. Such matter fields are termed 
as dark energy. Although we do not know at present what the dark energy is, there are many possible candidates 
proposed. The first and simplest model is a cosmological constant. Phenomenologically, the perfect fiuid model with 
■^ ■ an equation of state p = wp is often adopted from a data-analysis point of view, where w might be constant or 
^^ \ time-dependent. On the other hand, the simplest model for varying dark energy from a physical point of view is again 
OO . a scalar field with flat potential or possibly some other dynamical fields with appropriate potential 0, [j| . We here call 
such scalar field models for dark energy quintessence. There arc many variants of these varying dark energy models. 
If we restrict ourselves to the evolution of the homogeneous and isotropic universe, the perfect fluid and the 
quintessence models of dark energy play basically the same role with equivalent model functions, which is the equation 
of state in the former and the potential in the latter. However, once we turn our attention to inhomogeneities and/or 
J^ ' anisotropy, these two classes of models may show significant differences. Moreover, the response to inhomogeneous 
perturbations may distinguish the models which are degenerate in the homogeneous and isotropic evolution. Hence, 
it is indispensable to study inhomogeneities to distinguish the dark matter models. Our main interest in this paper 
is in the interaction between dark energy and black holes. 

The problem of mass accretion onto black holes in an expanding universe was raised by Zel'dovich and Novikov [5|, 
where they argued that the black-hole mass could increase self-similarly in proportion to the cosmological time. 
Although their argument was based on Newtonian gravity, self-similar solutions also arise in general relativity due to 
the scale-free nature of the Einstein field equation. Self-similar solutions are essentially characterized by functions of 
z = r /t and can describe inhomogeneous dynamics. They are also physically relevant because they may describe the 
asymptotic behavior of more general solutions. This is called self- similarity hypothesis [6| and, in fact, this was shown 
to be the case in some spherically symmetric gravitational collapse [7| ■ Sec [8|, |9| for a recent review of self-similar 
solutions and self-similarity hypothesis. See also [lO[ for a review of self-similar solutions in a more general context. 
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As we have a static black-hole solution in the Minkowski background, which is static, it would be natural to expect 
that we may have a self-similar black hole in the power-law flat Friedmann background, which is self-similar. In the 
study of the growth of primordial black holes [lH , Carr and Hawking [I4I and subsequent authors [ij, [ij, [I^, U3, 
found that, if we consider a perfect fluid with the equation of state p = (7 — l)p for 1 < 7 < 2, there are no self-similar 
solutions which have a black-hole event horizon and are asymptotic to the proper Friedmann solution at large distance, 
but there are self-similar solutions with a black-hole event horizon which arc only asymptotic to the Friedmann universe 
with some remaining anomaly. It has been realized that all the latter solutions arc only asymptotic to the Friedmann 
solution with anomaly in solid angle, which are termed as asymptotically quasi- Friedmann solutions [171. llSj. 

This historical problem has been recently revived by the discovery of the currently accelerated expansion of our 
universe. The accretion of dark energy or phantom energy onto a Schwarzschild black hole [i^,[20l and a Schwarzschild- 
de Sitter black hole [23 has been discussed. The cosmological evolution partially taken into account, it was suggested 
that black holes may grow self-similarly due to the accretion of a scalar field with potential [21[ . When the cosmological 
evolution is fully taken into account, however, it was shown [l6[ that there is no self-similar black-hole solution which 
is asymptotic to the decelerated Friedmann universe for a massless scalar field and a scalar field with positive potential. 
On the other hand, it has been recently found [23| that there is a one-parameter family of self-similar solutions which 
have a black-hole event horizon and arc asymptotic to the proper accelerated Friedmann universe for a perfect fluid 
with p = (7— 1)(0 (0 < 7 < 2/3). This strongly suggests that black holes can significantly grow due to the constant-rate 
accretion of dark energy in an accelerating universe. However, it should be noted that this phcnomcnological perfect 
fluid model for dark energy is ill-behaved in small-scale physics [l^. For a scalar field with such a flat potential that 
accelerates the Friedmann universe, it is still an open problem whether there is a self-similar black- hole solution which 
is asymptotic to the Friedmann universe. To answer this question, it is necessary to understand the properties of 
asymptotically Friedmann self-similar solutions containing a scalar flcld with potential, and this is investigated in the 
present paper. In spite of the motivation for self-similar black holes in the universe, the result obtained here generally 
applies to any objects which evolve in a self-similar manner and are embedded into the Friedmann universe containing 
a scalar field with potential. 

This paper is organized as follows. In Sec. II, we present a general formulation for self-similar solutions containing a 
scalar field with potential. In Sec. HI, we rewrite the field equations for nonlinear perturbations from the Friedmann 
solution. In Sec. IV, we find two independent one-parameter families of asymptotic solutions which are asymptotic to 
the Friedmann universe in different ways. In Sec. V, we present the physical properties of these asymptotic solutions. 
In Sec. VI, we summarize the paper. We use the units, in which c = 1. 

II. SELF-SIMILAR SOLUTIONS WITH A SCALAR FIELD 

We consider a single scalar field ip with potential V{ip) as a matter field, whose stress-energy tensor is given by 

Tab = f.affi - gab i ■^Z'-P .c^'"' + V {ip) j . (2.1) 

As we will see later, this can accelerate the expansion of the Friedmann universe. We adopt general relativity as a 
theory of gravity. The Einstein equation for this system is given by 
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v>,a^,b - gab [ -^y^.c^p'" + V{ip) 



(2.2) 



where k = vSttG and the comma denotes the partial derivative. The equation of motion for the scalar field is given 

by 

dip 

where D denotes the d'Alembertian associated with gab- We consider a spherically symmetric spacetime, in which the 
line element is given by 

ds^ = -e2*(*''-)di2 + e2*(*''-)dr2 + R\t, r)dn\ (2.4) 

where dVP' = dO^ + sin^ Odcj)^ is the line element on the unit sphere and the domain of 9 and (j) arc Q < 9 < tt and 
< (/)< 27r. 

We assume that the spacetime is self-similar, which is defined by the existence of a vector field £^°' such that 

^igab = 2gab, (2.5) 



where C^ denotes the Lie derivative along ^''. This vector field ^° is called a homothetic Killing vector. If ^° is tilted 
to {d/dt)°-, nondimensional metric functions depend only on z = r/t [2J| , i.e., 



$ = $(z), * = *(z), R = rS{z). 
Then, the scalar field tp and its potential V{(p) are of the following form |25| : 



(2.6) 



if = — -lnr + /(z). 
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where Vq and A are constants. It should be noted that for simplicity we have assumed that r and t are positive. In 
fact, we can recover the results for the general case simply by replacing r, t and z with |r|, |i| and \z\, respectively. 

Since we are interested in self-similar solutions perturbed from the Friedmann universe at large distances, we can 
assume that the gradient of the scalar field is timelike near spacelike infinity. In such a case, we can take time slicing 
so that (p depends only on the time coordinate t, which we call the constant scalar field slicing. In fact, this coordinate 
system is equivalent to the comoving coordinates where only the diagonal components of the stress-energy tensor are 
nonvanishing. In the following we choose this slicing, so that we have 



f{z) = -Inz 
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(2.9) 



and 






(2.10) 



where ipo is a constant. If the scalar field is massless, i.e., Vq = 0, we can simply delete (po because only the gradient 
of (f appears in the action. If the scalar field has a potential, we can renormalize the constant ipo in Eq. (j2.10p into 
the factor Vq in the scalar field potential by replacing Vq with Vq such that 
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Therefore, we set (po = in the following. 
In this coordinate system, tt, tr, rr, and 



components of the Einstein equation, respectively, yield 
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where the prime denotes the ordinary derivative with respect to In z and 



K = ze*-* 



(2.16) 



is the relative velocity between the constant z surface to the constant r surface. The equation of motion for the scalar 
field becomes 



^V? ($'-*'- 2C - 1 
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Four of the five equations (|2.12p - (|2.15p and (|2.17p are independent 
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(2.17) 



We derive the following two relations for later use. Adding Eq. (I2.12p to Eq. (I2.14p and using Eq. (|2.13p . we get 

T/2^'-$' = - — V;2 ^2.18) 

A 

Subtracting Eq. ((TT^ from Eq. ([^1^ and using Eq. (j^J^ . we get 



III. NONLINEAR PERTURBATION FROM THE FRIEDMANN SOLUTION 

A. The flat Friedmann solution in self-similar coordinates 

The flat Friedmann spacetime is given by the following line element: 

ds^ ^^dt^ + a{t)^{df^+f^dn^). (3.1) 

The scale factor a{t) satisfies the Friedmann equation 

^)"4»^(i.' + n.)). (3-2) 

where the dot denotes the derivative with respect to i. (p = ip{t) satisfies the equation of motion 

>fi + 3Hip+-}^ = 0. (3.3) 
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For these equations, we have a power-law solution 



a = aot", (3.4) 

ip = -^Int, (3.5) 



where Oq is a constant and a and Vq are given by 
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This is obviously compatible with Eq. (|2.10p . A massless scalar field formally corresponds to A^ = 6, where a ~ 1/3 
and the cosmic expansion is decelerated. For a nontrivial potential, if < A^ < 2, then the cosmic expansion is 
accelerated, while, if A^ > 2, then the cosmic expansion is decelerated. For A^ > 6, the potential becomes negative. 
If a 7^ 1 , relating r to f through 

aof=^, (3.8) 

we can rewrite the flat Friedmann solution in the standard form for self-similar spacetimes, where 

e* = l, e* = z-", S^-^--z-''. (3.9) 

|l-a| 

So, the power-law flat Friedmann solution is self-similar. If a = 1, the flat Friedmann solution is still self-similar but 
the homothetic Killing vector is parallel to {d/dtY (see e.g. [26[). This case needs a special treatment and we do not 
consider this case in the present paper. It should be emphasized that from Eq. p.8p f ^ and f — > c» correspond 
to r — > and r — > c», respectively, for < a < 1, while this is reversed for a > 1. When we study spatial infinity in 
general case, we should take the limit z^~" -^ oo for fixed t [27| . 



B. Field equations for nonlinear perturbation 

Since we are interested in self-similar solutions which are asymptotic to the Friedmann solution, we write general 
spherically symmetric self-similar solutions in the following form: 

e* = e^W, e* = z-"e^(^\ ^ = ^^z-"e^("\ <^ = A inf + Z?(z). (3.10) 

|1 — aj kX 

As for the gradient of the scalar field, we get 
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(3.11) 



The above equation implies that, if A and B arc finite and D' is sufficiently small, we can choose the constant scalar 
field slicing, where D is a constant Dq. This is the case where the solution is asymptotic to the flat Friedmann 
solution. Hereafter we take the constant scalar field slicing. We can set Dq = 0. Then wc have 

e* = e^(^), e* = z-"e^(^), S = -^z-''e^'^^\ ^^\lnt. (3.12) 

|1 — a\ kX 

Substituting the above, we derive the following set of ordinary differential equations for A, B and C": 

- aA' + A'C + (1 - a)B' + B'C - C" - (1 - a)C' - C"^ = (3.13) 

from Eq. ((TTO)) . 

A' -B' - 2C" - (3a - 1) {e"^^ - 1) = (3.14) 

from Eq. (PTZ)) . 

V^B' -A'^0 (3.15) 

from Eq. (j2.18|) . and the constraint equation 

(142 _ i)(c"2 _ 2aC" - aA') - 2C' + [a - lf{e^^-^^ - 1) = (3.16) 

from Eqs. ^JQ, (PH)) and ([XT5)) . where V^ is given by 

y^2 ^ ^2-2ag2B-2A (3^^7) 

from Eq. (|TTC)) . 

IV. ASYMPTOTICALLY FRIEDMANN SOLUTIONS 

Although we are most interested in asymptotically proper Friedmann solutions, wc also study more general solutions 
to which the asymptotic scheme applies. So we only require that all A, B and C have finite limit values, i.e., 

A^Ao,B^Bo,C^Co (4.1) 

at spatial infinity, i.e., as z^~" — > 00. Then, A' , B' , and C" tend to vanish from rHospital's rule. Hereafter, we always 
choose Aq = Bq = hy rescaling the coordinates t and r, whereas Cq may not vanish. 

It is not so trivial how such asymptotic solutions are expanded around z^~" — > cx). First, we note that only V~'^ is 
explicitly of higher order in Eqs. (|3.13|) - (|3.15|) in the present limit. Equation (|3.15p then implies that B' is always of 
higher order than A' . If wc linearize Eqs. (|3.13p and (|3.14|) . we get 

- aA' - C" - (1 - a)C' = 0, (4.2) 

and 

A' - 2C' - 2(3a - 1)A = 0, (4.3) 



respectively. Then, eliminating C" and C", we get 

A" - (5a - 3)^' + 2(a - l)(3a - 1)^ = 0. (4.4) 

A general solution of the above equation is given by the linear combination of the following two independent solutions: 

A = ^^"-1 (4.5) 

and 

A = z^°'-^. (4.6) 

The second solution is always valid in the limit z^~" — > cx), while the first is valid only for {ia — 1)/(1 — a) < 0. Since 
a is positive, the first solution is valid only for < a < 1/3 or a > 1. 

A. Asymptotically proper Friedmann solutions 

To get the full form of the first solution (|4.5p . we use Eqs. p.l3p - p.l5p and the result is the following: 

A = Axz^"-^, (4.7) 

B = Biz^"-^, (4.8) 

C = Co + Ciz^"-^ (4.9) 

in linear order, where 

„ 3a — 1 . , , 

Bi = -A,, 4.10 

5a — 3 

C\ = -\a^. (4.11) 

It turns out that we need higher order terms to see whether we can have nontrivial solutions and whether Co 
vanishes or not in Eq. p.l6p . It is cumbersome but straightforward to get higher order terms from Eq. (|3.13[l - (|3.15p . 
The result is 

A = Aiz^"-^ + A2z'^"^^ + Asz^"-^ + ■ ■ ■ , (4.12) 

B = Biz^"-^, (4.13) 

C = Co + Ciz^"-! + C2z^""^ + Cgz^""^ + • • • , (4.14) 



where the coefficients are all parametrized by Ai as follows: 



^2 = -^A,, (4.15) 

a — 1 

4a2 - 3a + 1 , , , 

^^ - 2(5a^3)(l-a) -^^' ^'-''^ 

11a — 1 .n , , 

A3 = -^^Al (4.17) 

Cs = -^Al (4.18) 

For a > 1, we consider the limit z — *■ 0. In this case, the lowest order of Eq. (j3.16p . which is of order z°, just 
yields Co = 0. Therefore, we get self-similar solutions which are asymptotic to the proper Friedmann solution with 
vanishing Cq. These solutions are termed as asymptotically proper Friedmann solutions. 

For < a < 1/3, the situation is more complicated. In this case, we consider the limit z -^ oo. Then, the terms 
of order z^"-^ in Eqs. (|iT^ and (|iT^ get higher than those of order z^"-^. Substituting Eqs. ([JT^ - fiH)) into 
Eq. (|3.16p . and using Eqs. (|4.17p and (|4.18p . we can see that the terms of orders 2^+" and 2^" all cancel out. Also 
in this case, the nontrivial lowest order, which is of order z°, yields Co = 0. So, these self-similar solutions are 
asymptotically proper Friedmann solutions. 

For a = 1/3, from the linear order analysis, we find that A = Ai and C = Ci are constants, while B vanishes from 
Eq. (|4.10p . It should be noted that we may have higher order terms. As we have set Aq = 0, we can set Ai = by 



rescaling the time coordinate. Then, if we have higher order terms, they must satisfy Eq. (|4.4p and this again yields 
A = const and A = z^"~^. The latter case must be included into the next case. Hence, we can concentrate on the 
solution where A = 0, B ^ and C — Ci. In this case, we can show Ci = from Eq. (|3.16p . Therefore, the solution 
coincides with the exact Friedmann solution. 

In summary, there is a one-parameter family of asymptotically proper Friedmann self-similar solutions for < a < 
1/3 or 1 < a. There is no nontrivial asymptotically proper Friedmann self-similar solution for 1/3 < a < 1. 

B. Asymptotically quasi-Friedmann solutions 

Up to the nontrivial lowest order, the second solution (|4.6p is given by 

A = Aiz^''-^ 
B = Biz'^"-'^, 
C = Co + Ciz2"-2, (4.19) 

where 

Bi = Ai, 

Ci = T^M. (4.20) 

\ — a 

From the lowest order of Eq. (|3.16p . which is of order z°, we get 

Higher order terms are expanded in terms of integer powers of z^"^" and the coefficients are written by integer power 
of Ai. Hence, if Cq = 0, we have Ai = Bi = Ci = and the solution becomes trivial. Only if Cq 7^ 0, we have 
a nontrivial solution. We term these nontrivial solutions as asymptotically quasi-Friedmann solutions. There is a 
one-parameter family of such solutions for < a < 1 or 1 < a. 

V. PHYSICAL PROPERTIES OF THE SOLUTIONS 

We have shown that there are two types of asymptotically Friedmann solutions with trivial and nontrivial asymptotic 
values for C . The first is asymptotically proper Friedmann and the second is asymptotically quasi-Friedmann. In this 
section we see their physical properties. 

A. Solid angle anomaly 

We have the following asymptotic form of the metric near spatial infinity: 

-2a 

-dt^ + alt'"idf^+e'''°r'dn^), (5.2) 



ds^ - -dt^ + z-^'^di'^ + y^ ^e^^'or^dn^ (5.1) 



where Cq = and Cq ^ hold for asymptotically proper Friedmann and quasi-Friedmann solutions, respectively. 
If we consider a two-sphere on which t = const and r = const, its area is given by 47ra§i^"f^e^'^'', while the proper 
length of the radius on the constant t hypersurface is equal to aof^f. So the ratio of the area to the squared radius 
is not 47r but 47re^'^°. Therefore, there is a surplus in the solid angle for Co > and a deficit for Co < 0. Only for 
Co = 0, we have no anomaly in the solid angle. 

We can see this metric in another way. When we consider the 9 = tt/2 section, we get the line element 

ds^ = -dt^ +a^t^°'{df^ +e'^^'>f'^d(l)^) (5.3) 

= -dt^ + a^t^°'{df^ + r^d4)^), (5.4) 



where (/) = e "(/) and hence 

< < 2iTe^°. (5.5) 

Although the above Hne element is the same as that for the = 7r/2 section of the Friedmann solution, the domain 
of the azimuthal angle is anomalous. In fact, there is a surplus in the azimuthal angle for Cq > and a deficit for 
Co<0. 

This kind of anomaly in the solid angle is already discussed in the context of static global monopoles and termed 
as solid angle deficit for Cq < [2^. Hence, we can say that the asymptotically quasi- Friedmann solutions are with 
solid angle surplus or deficit, while the asymptotically proper Friedmann solutions are not. It should be noted that 
despite the apparent similarity with conical singularities in cylindrically symmetric spacetimes, the spacetime with 
the solid angle anomaly is not flat even locally. 

B. Density perturbation 

It is also interesting to get insight into the difference of the two classes of asymptotic solutions in the density field 
at spatial infinity on the constant t hypersurface. The energy density p observed by a comoving observer is given by 



p = naUbT'^'' = i 
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(5.6) 



where n° is a unit vector normal to the constant scalar field hypersurface. Hence, the background Friedmann density 
Pb, the density perturbation Sp, and the density contrast Ap are, respectively, given by 

Pb = 7^2^^' (5-7) 

5p = p{t,r) - p^{t,r) ^ --^^^A, (5.8) 



Pb 3 



Ap ^ ^ « -^A, (5.9) 



where the suffix b denotes quantities for the background Friedmann solution and the weak equality "sa" denotes that 
the ratio of both sides approaches unity in the relevant limit. The asymptotic form of the physical areal radius R is 
given by 

R^rS^ , ^ rz-^e^" (5.10) 

|l-a| 

for both cases. So, the fall-off of the density perturbation in terms of the physical areal radius is given by 

3q-1 

and 

5p ex -Aii"2^2a-2 ^ -AiR-^ (5.12) 

for asymptotically proper Friedmann and quasi-Friedmann solutions, respectively. Therefore, the density perturbation 
rapidly falls off for asymptotically proper Friedmann solutions for the accelerated case a > 1. It falls off as R^^ for 
asymptotically quasi-Friedmann solutions. For asymptotically proper Friedmann solutions with < a < 1/3, where 

the potential is negative, the fall-off is as slow as R^ ^ ~ on the constant t slice, which is much slower than 

for asymptotically quasi-Friedmann solutions. 

C. Mass perturbation 

In spherically symmetric spacetimes, the Misner-Sharp mass m is known to be a well-behaved quasilocal mass 
defined as [29l.[3(j| 

^ ^ 1 + e-2*i?2 _ e-2*i?2^. (5.13) 



In the present formulation, this can be rewritten as 



y = 1 + 7r^7^^'"'"e2C^"'-'(-" + C")' - (l + Y^) ' e^^-^^, (5.14) 



where M is the nondimensional mass defined by 



M=—. (5.15) 

r 



For the flat Friedmann solution, this quantity becomes 



z^'^'-^a^(^\ . (5.16) 



2"A a^ 2~2a _ „2 f Rh 



The perturbation for this quantity 



fi/b (l-«) 



M^^hr (-)-hr W (5-17) 



it/ \ R J yi?/j^ 



is given by 



'^ f'^) " -Jl^"^'"'^" ^ -"'' - "l*'""'^^^'""^^i (t) ^'"""^^^''"^ (5.18) 



for asymptotically proper Friedmann solutions and 

for asymptotically quasi-Friedmann solutions. Hence, the ratio of the perturbation to the background value 



is given by 



and 



A«-iAiz3"-i = _^^ Ai(i^) (5.21) 

a a \ t J 



A « e^'^" - 1 (5.22) 

for asymptotically proper Friedmann and quasi-Friedmann solutions, respectively. 

It should be noted that since m and R are nonlinearly perturbed for asymptotically quasi-Friedmann solutions, this 
ratio depends on whether we compare the perturbation at the same r/t or at the same R/t. Noting that S is also 
perturbed, the mass perturbation SM{z) = M{z) — M-t,{z) is directly given by 

Am = ^ = (1 + A)e^ - 1. (5.23) 

Mb 



For asymptotically proper Friedmann solutions, Aj\/ is calculated as 



(3Q-l)/(l-a) 



Am « -^^1^'""' « -^1 (t^) • (5-24) 

Hence, the relative mass perturbation Aj\/ tends to vanish at spatial infinity for both a > 1 and < a < 1/3. This 
also implies that the mass perturbation Sm itself tends to vanish for a > 1 but diverge for < a < 1/3 as it is 
proportional to (i?b/i)^^*-^~"-'i- For asymptotically quasi-Friedmann solutions. Am is calculated as 

Am ~ e^'^" - 1. (5.25) 
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Hence, the relative mass perturbation tends to be constant. This is directly related to the solid angle anomaly. If 
A]^J is positive (negative), there is a solid angle surplus (deficit). This situation is apparently opposite to the case of 
global monopoles, where the positive (negative) mass density implies deficit (surplus) in the solid angle. This is due 
to the fact that the Misner-Sharp mass is dominated by the first and third terms on the right-hand side of Eq. (|5.13p 
for the static configuration, while it is by the second term, i.e., the kinematic term, for the flat Friedmann solution. 
The mass perturbation 6m itself diverges as {R/t)^t for both < a < 1 and 1 < a. 

So, in order to have an asymptotically proper Friedmann solution from an accelerated Friedmann universe, for which 
the potential is positive, we only need to perturb a finite amount of mass and the mass perturbation is compensated at 
spatial infinity. In contrast, in order to have an asymptotically quasi-Friedmann solution from the Friedmann solution, 
we need to perturb an infinite amount of mass and the mass perturbation remains at spatial infinity. This is also the 
case so as to have an asymptotically proper Friedmann solution from the decelerated Friedmann solution, for which 
the potential is negative. This suggests that asymptotically proper accelerated Friedmann solutions are physically 
acceptable as nonlinearly perturbed solutions from the Friedmann solution by some classical mechanism. This also 
suggests that any classical perturbation mechanism will not perturb a Friedmann universe to a quasi-Friedmann 
universe. Only through quantum fiuctuations, it might be possible to have a quasi-Friedmann solution because an 
infinite amount of perturbed mass must extend in scales much larger than the Hubble horizon at any epoch. On the 
other hand, whether the perturbed mass is compensated or remains or even diverges at spatial infinity, the present 
perturbation scheme is still completely applicable for these asymptotic solutions. 

D. Comparison ^vith a perfect fluid with p = (7 — l)p 

For the Friedmann solution, a scalar field with exponential potential and a perfect fluid with p = (7 — l)p play a 
completely equivalent role. In the spatially flat case, they are related with the following relation: 

«=A.^. (5.26) 

So, the accelerated expansion is possible if < A^ < 2 for the scalar field and if < 7 < 2/3 for the perfect fluid. 
However, once we admit perturbations from a uniform distribution, the two systems get very different. 

For example, a scalar wave propagates at the speed of light in the short wave length limit in the scalar field system 
even in the presence of potential. In contrast, in the perfect fiuid system with the equation of state p = (7 — l)p, a 
sound wave propagates at the sound speed V7 — 1 for 1 < 7 < 2 and, in fact, there is no sound wave but instability 
in the short wave length limit for < 7 < 1 [18l |. 

Also in the perfect fiuid system with p = (7 — l)p, there are two independent one-parameter families of solutions 
which are asymptotic to the Friedmann solution |17l . Il8j . One is asymptotically proper Friedmann solutions at spatial 
infinity and the other is asymptotically quasi-Friedmann solutions at spatial infinity. The latter is valid for both the 
accelerating (0 < 7 < 2/3 or a > 1) and decelerating (2/3 < 7 < 2 or 1/3 < a < 1) cases, while the former is only 
valid for the accelerating case. Hence, the situation is exactly parallel to that in the scalar field case. This is a very 
unexpected result because we do admit inhomogeneity when we consider asymptotic solutions. 

In the perfect fiuid analysis, the strongly decelerated case, < a < 1/3 or 7 > 2, has not been analyzed because 
causality is violated in such a model. In the present analysis, on the other hand, since the scalar field with negative 
potential is acceptable from a causal point of view, we have included this case and found interesting features that both 
asymptotically proper and quasi-Friedmann solutions exist and that the asymptotically proper Friedmann solutions 
are very different from those for the accelerated case. 

For a perfect fluid with 1 < 7 < 2, there is no self-similar solution which has a black-hole event horizon and is 
asymptotic to the proper Friedmann solution at spatial inflnity However, for a perfect fluid with < 7 < 2/3, there is 
a one-parameter family of asymptotically proper Friedmann solutions. In fact, the numerical integration has revealed 
that there is a one-parameter family of self-similar solutions among them which contain a black- hole event horizon [23] . 
To implement the numerical integration in that case, it is highly advantageous that the system of ordinary differential 
equations has no critical surface because there is no propagation of sound wave. Also in the scalar field case, one 
might guess that the existence of asymptotically proper Friedmann solutions suggests the existence of self-similar 
black-hole solutions belonging to this class. However, the scalar field system has a critical surface coinciding with 
a similarity horizon, where Vz = 1. This makes the problem complicated because this could possibly increase the 
number of self-similar solutions drastically as such a critical surface may admit weak discontinuity. In this connection, 
we should also note that because of the critical surface, the power-law flat Friedmann solution containing a scalar field 
with potential is unstable for 4 < A^ < 6 against weak discontinuity, i.e., the kink mode at a particle horizon [3l|. It 
is however stable for < A^ < 4, marginally stable for A-^ = 4. This kink instability might be related to the physical 
relevance of self-similar solutions. 
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VI. SUMMARY 

We have considered self-similar nonlinear perturbation from the Friedmann solution and investigated the asymptotic 
properties of spherically symmetric self-similar solutions containing a scalar field with potential and approaching the 
flat Friedmann solution at spatial infinity. The potential is restricted from self-similarity to be exponential with the 
steepness parameter A. This is motivated by the fact that a scalar field with sufficiently flat (0 < A^ < 2) potential 
enables the universe to expand with acceleration and hence acts as quintessence. 

If the potential is so flat, i.e., < A^ < 2 that the Friedmann universe expands with acceleration, we have found 
that there is a one-parameter family of self-similar solutions which are asymptotic to the proper Friedmann solution 
at spatial infinity. Furthermore, we have found that there is also a one-parameter family of self-similar solutions 
which are asymptotic to the Friedmann solution but with some anomaly in solid angle. Such solutions are called 
asymptotically quasi- Friedmann solutions. 

If the potential is steep, i.e. A^ > 2, we have the Friedmann universe decelerated. Even in such a potential, we have 
found a one-parameter family of self-similar solutions which are asymptotically quasi-Friedmann solutions. However, 
we have also shown that there is no nontrivial asymptotically proper Friedmann self-similar solution in this case as 
long as the potential is positive. We should note that it was already shown that there is no self-similar solution which 
contains a black-hole event horizon and is asymptotically proper Friedmann or quasi-Friedmann for a scalar field with 
positive potential inducing the decelerating expansion [l6| . 

Our analysis includes the case of a massless scalar field, where the flat Friedmann universe is decelerated. In this 
case, we have found that there is a one-parameter family of asymptotically quasi-Friedmann self-similar solutions, 
while there is no nontrivial asymptotically proper Friedmann self-similar solution. We should also note that it was 
already shown that there is no self-similar solution which contains a black-hole event horizon and is asymptotically 
proper Friedmann or quasi-Friedmann for a massless scalar field [16[ . 

Our analysis also includes the case where the potential is negative. In such a case, the Friedmann universe is 
strongly decelerated. We have found that there are both one-parameter families of asymptotically quasi-Friedmann 
self-similar solutions and asymptotically proper Friedmann self-similar solutions. The latter is very different in density 
and mass perturbations from that for the positive potential. 

We have shown that the perturbed mass is finite for asymptotically proper Friedmann solutions as long as the 
potential is positive. In contrast, it is infinite for asymptotically quasi-Friedmann solutions. This suggests that 
asymptotically proper Friedmann solutions are physically more acceptable as solutions perturbed from the Friedmann 
universe through some causal mechanism than asymptotically quasi-Friedmann solutions. Although asymptotically 
proper Friedmann solutions are possible even if the potential is negative, the perturbed mass is infinite there. 

Although we have found the above interesting properties of self-similar solutions containing a scalar field with 
potential, it is still an open question whether there is a self-similar black-hole solution which is asymptotically proper 
or quasi-Friedmann. We need possibly a numerical analysis based on the present asymptotic analysis, as it has 
revealed the existence of self-similar black-hole solutions for a perfect fluid with the equation of state p = (7 — l)p 
(0 < 7 < 2/3) [l^, [2^. Although spherically symmetric self-similar solutions with scalar fields have been also 
investigated in a dynamical systems approach |32l . ISSl . |3J| , no definite answer to the existence of black-hole solutions 
has been reported yet. It is an important future work to answer whether there is a self-similar black- hole solution in 
quintessential cosmology and, if it exists, to study the physical properties of such a black-hole solution. 
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